Are f(R) dark energy models cosmologically viable ? 
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All f{R) modified gravity theories are conformally identical to models of quintessence in which 
matter is coupled to dark energy with a strong coupling. This coupling induces a cosmological 
evolution radically different from standard cosmology. We find that in all f{R) theories that behave 
as a power of R at large or small R (which include most of those proposed so far in the literature) the 
scale factor during the matter phase grows as t 1 ^ 2 instead of the standard law t 2 ^ 3 . This behaviour is 
grossly inconsistent with cosmological observations (e.g. WMAP), thereby ruling out these models 
even if they pass the supernovae test and can escape the local gravity constraints. 
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The late time accelerated cosmic expansion is a ma- 
jor challenge to cosmology It can be due to an ex- 
otic component with sufficiently negative pressure, Dark 
Energy (DE) , or alternatively to a modification of grav- 
ity, no longer described by General Relativity. Examples 
of such modified gravity DE models are theories where 
the Ricci scalar R in the Lagrangian is replaced by some 
function f(R), e.g. inverse powers R~ n 0, 0]- Although 
these models exhibit a natural acceleration mechanism, 
criticisms emphasized their inability to pass solar sys- 
tem constraints Q. Indeed, f{R) theories correspond 
to scalar-tensor gravity with vanishing Brans-Dicke pa- 
rameter wbd 0- However one could in principle build 
models with a very short interaction range (e.g. adding 
a R 2 term 0, 0|) or assume decoupling of the baryons 
from modified gravity. Since these models could pass 
local gravity constraints, it is important to assess their 
cosmological viability: this is the aim of this Letter. We 



will consider models of the form f{R) = R — jj, 



2(71+1) 



/R n 



for all fx,n such that df/dR > 0; for all these models, 
the scale factor a(t) expands as i 1 / 2 instead of the con- 
ventional £ 2 / 3 behaviour during the matter phase that 
precedes the final accelerated stage (in contrast with in- 
flationary models like Starobinsky's R 2 one This 
would lead to inconsistencies with the observed distance 
to the cosmic microwave background (CMB), the large 
scale structure (LSS) formation, and the age of the Uni- 
verse. This crucial fact appears to have been overlooked 
so far. 

Consider the general action in the Jordan frame (JF) 



S = 



2k 



;f(R) + £r, 



(1) 



where k 2 = 8nG (G is the gravitational constant). For 
a flat Friedmann-Robertson- Walker (FRW) metric the 
equations are given by 



3FH 2 = (RF - f)/2 - 3HF - 
2FH = —F 



HF - k 2 ( Pn 



Pr, 



Pm 1 

,), (2) 

where F = df/dR, H = a/a, and p m and p m represent 
the energy density and the pressure of a perfect fluid, 



obeying the standard conservation equation. These equa- 
tions coincide with a scalar-tensor Brans-Dicke theory 
with a potential and vanishing wbd @, [13] • 

Under the conformal transformation g^ = e^g^, 
2uj = logF, one obtains the Einstein frame (EF)action: 



Se = 



d 4 x 



m 

2k 2 



■(V0) 2 -V(<f>)+C m (<f>) 



,(3) 



where (j> = VGu/k and V = sign(F){RF - !)/2k 2 F 2 
(all tilded quantities are in EF). The conformal trans- 
formation is singular for F = 0, so we will consider only 
positive-definite forms of F. Quantities in the two frames 
are related as follows 



-4cu 



p m e 



-4a) 



at 



J dt, 



J a.(4) 



Although we will work mainly in the EF, we checked 
all numerical and analytical results directly in the JF as 
well. In EF the field <f> and the fluid satisfy the standard 
gravitational and conservation equations: 

$ + 3Hi + Vj = y/WKP(fim-3Pm): ( 5 ) 
p m + 3H (p m + Prn) = - \f2~j3 Kf3(p(p m - 3p m ), (6) 



where the coupling (3 is given by 
= 1/2, 



(7) 



regardless of the form of f(R). Then the strength of the 
coupling between the field and the fluid is uniquely de- 
termined in all f(R) gravity theories. A dimensionless 
strength of order unity means that matter feels an addi- 
tional scalar force as strong as gravity itself. Note that j3 
is related to wbd via the relation (3 = [3/4(2wbd + 3)] 1 / 2 . 
The dynamics of the system depends upon the form of 
the potential V{4>), i.e., the choice of f(R). For theories 
in which f(R) = — p 2 ^ n+1 ^ R~ n (n ^ —1, 0, negative n are 
also included), the potential in EF is a pure exponential 



V{4>) = Acxp {-\K<j>) 



(8) 
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Figure 1: Evolution of the fractional energy densities (<f>, mat- 
ter, radiation) in EF for the model f(R) = R — /J- 4 /R (top 
panel). Overimposed as dotted lines the evolution of for 
n = 4 and 10. Notice the constant value £7^ ~ 1/9 in the 
itjMDE phase between radiation and DE domination. In the 
bottom panel we plot the evolution of the observed EOS iode 
of DE in JF and the effective EOS in both EF and JF (n = 1). 



are constant and given by (ill . fl2| fi^ = w e g = 4/3 2 /9 . 
Then we have fi^ = uW = 1/9 in /(i?) gravity theories, 
regardless of the form of f(R). Therefore, contrary to 
standard cosmology, in coupled models DE is not neg- 
ligible in the past (until the radiation era). In contrast 
to the accelerated attractor, the 0MDE occurs when the 
coupling term in the r.h.s. of Eq. |(5]) dominates over the 
potential term, as it can be explicitly shown. This aspect 
is crucial for the present work since it implies that the 
0MDE exists independently of the form of f(R). In this 
regime the scale factor behaves as a ~ £ 3 / 5 : in the JF 
this becomes a ~ t 1 / 2 instead of the usual t 2 / 3 behav- 
ior of the MDE. This is clearly a strong deviation from 
standard cosmology and, as one can expect, is ruled out 
by observations, as illustrated below. Notice that the JF 
evolution in this phase corresponds to R = as during 
the radiation epoch but, just as in that case, there is no 
singularity in an inverse power-law theory because this 
behavior is not exact (see below). In the language of 
dynamical systems, the 0MDE is a saddle point. 

To analyse numerically the full system (including the 
radiation energy density p ra( j which obeys the standard 
conservation equation in both frames), we introduce the 
following quantities: 



x\ 



k4>' 

7f : 



x 2 



K 

H 



X3 



K 

H 



Prad 



(9) 



where A 



y/E n+2 



and A 



M 2 



• The condition 



3 n+l " 2n 2 n\n\ 1 

F > implies A > except for — 1 < n < 0: in this case 
since the potential becomes negative we analyse directly 
the JF. In EF the R~ n model corresponds to a coupled 
DE scenario studied in Refs. [ill . [l2| with the coupling 
(O. We first discuss the main properties of this expo- 
nential potential and then extend them to the general 
case. 

As shown in pjj for all values of n outside (—1,0) 
the system has one and only one global attractor solu- 
tion, a scalar-field dominated solution with an energy 
fraction 0^ = 1. This solution appears when the poten- 
tial term in Eq. ((5J dominates over the coupling term 
on the r.h.s., and is therefore independent of the cou- 
pling. On this attractor the scale factor evolves as 
a ~ p/[3(i+t«cff)] w here the effective equation of state 

(EOS) is w cS = -1 + A 2 /3 = -1 + l^XS ■ This can be 
identified with the acceleration today if /x ~ Ho ■ 

Beside the final attractor, a coupled field with an expo- 
nential potential has also another solution in which mat- 
ter and field scale in the same way with time and, conse- 
quently, their density fractions are constants. This epoch 
has been denoted as (^-matter-dominated era (</>MDE) 
[ll[. As we will show in a moment, the </>MDE plays 
a central role in this work. This epoch occurs just after 
the radiation era and replaces the usual MDE. During the 
0MDE the energy fraction and the effective EOS w c s 



where a prime denotes the derivative with respect to TV = 
log (a). The energy fractions of the field (j> and of matter 
are given by = x\ + x\ and f2 TO = 1 — x\ — x\ — x\ , 
respectively. The effective EOS and the field EOS are 
w c ff — x\ — x\ + x§/3 and = {x\ — X2)/(x\ + x\), re- 
spectively. The complete system has been already stud- 
ied in Ref. [n| and we will not repeat it here. The (TjMDE 
corresponds to the fixed points {x\,X2,Xs) — (1/3,0,0) 
with Vttf, = w c ff = 1/9. After this, the universe falls on 
the final attractor, which is the accelerated fixed point 
(xi,X2,X3) = (A/ \/6, y/l — A 2 /6, 0) with 0^ = 1 and 
w cS = -1 + A 2 /3. 

To return to the JF we can simply apply the transfor- 
mation law (@]). In the regime where radiation is negligi- 
ble (X3 ~ 0) we obtain the following effective EOS 



2n 



3(1 + n) (l-xi) 2 



together with the relation 



K 2 pr, 



3H 2 F (l-xi) 2 ' 



(10) 



(11) 



For the accelerated attractor we have w; e ff = — 1 + 
3(i+n){i+27t) (this relation was originally found in the 
context of inflation in [ll|), which gives w g = —2/3 for 
n = 1. The 0MDE corresponds to X2 = and therefore 
w c ff = 1/3 for any n, giving a cx T 1 / 2 . From Eq. ifTTj) one 
has n m = 2 and Or = (RF - f - 6HF)/6FH 2 = -1 



3 



in the 0MDE [see Eq. J5J]. Since does not need to 
be positive definite, Q m can be larger than unity without 
any inconsistency. 

Notice also that w^, differs from the quantity wde used 
to annalyse SN data and defined through the equation 
H 2 = #2[fi^ a - 3 + (1 - f2^)a" 3 exp(-3 / daiu DE /a)]. 
We will also computte »de below. 

Most modifications of gravity suggested in the liter- 
ature consider terms in addition to the usual Einstein- 
Hilbert Lagrangian. For instance, several authors have 
studied the following DE model 0, S| : 

f{R)=R-H 2(n+1) /R n . (12) 
In this case the potential in EF is given by 

V{(j>) = Ae-^P**^"* -1)^t, (13) 

which vanishes at <fi = and has a maximum at n(f) = 
2(n + l)/{n + 2). In the limit (j) — > oo it behaves as 
V(<p) oc exp(— Xnif). For negative n the pure exponen- 
tial approximation is always good during the past cosmic 
history if the higher-curvature term is responsible for the 
present acceleration, because we are always in large R 
limit. For positive n the potential (fl3|) differs from the 
pure exponential one in the limit R ^> /i 2 and one might 
expect that for these values the evolution goes on as in 
the standard case and, in particular, the </>MDE disap- 
pears. However, we show now by building an explicit 
solution that in reality this does not happen. 

Let us focus on the R — ^ 4 / R model taken for sim- 
plicity without radiation. During the 0MDE one can 
approximate the scale factor in JF as a(t) = (t/U) 1 / 2 + 
e(i)(i/*i) 9/4 where U is an initial time at the beginning of 
the </>MDE. This solution is valid at first order in e pro- 
vided e = (n 2 /U4Hf)/((p$/3H?) - ^H/H") 1 / 2 , which 
is indeed small for /i of order Hq as present accelera- 
tion requires. After some time, the correction gets larger 
than the zero-th order term itself and the 0MDE is fol- 
lowed by a phase of accelerated expansion. Then the 
beginning of the late-time acceleration is quantified by 

the condition t « (lUHf^J (p$ /3F?) / £t 2 ) 4/7 <<• A sim- 
ilar argument applies for any n < — 1, n > —3/4 with a 
correction growing as t 5 / 2 - 1 / 2 (™+ 1 ). Since R is of order 
/i 2 at the beginning of the </>MDE, the term /j, 4 /R domi- 
nates over R after the radiation era. This applies for any 
/x, no matter how small. In other words, the limit — ► 
of a fourth-order theory does not reduce to second-order 
general relativity if, at the same time, one imposes the 
conditions of acceleration today (cfr. [1J|). 

For larger [i (>> Hq) the </>MDE can be shortened or 
bypassed from the above condition, but then DE domi- 
nates soon without a matter dominated epoch. Thus we 
have only two cases: either (i) the </>MDE exists, or (ii) 
a rapid transition from the radiation era to the acceler- 
ating stage (without 0MDE) takes place. In summary, 
the system never behaves as in a standard cosmological 
scenario except during radiation (during which matter 



and field play no role in the expansion rate). In other 
words, whenever matter is dynamically important, DE is 
also important as a consequence of the coupling. We now 
confirm all this by a direct numerical integration. 

In Fig. Q] we plot the evolution of fi^, Odm and f2 ra d 
and the equation of state in EF for n = 1, 4 and 10 and 
/i « H . The present values of the radiation and matter 
density fractions are chosen to match the observations in 
JF. As expected, the system enters the 0MDE after the 
radiation era and finally falls on the accelerated attrac- 
tor. We ran our numerical code for other positive and 
negative values of n (from n = — 10 to n = 10, limiting 
to the accelerating cases) and found similar cosmological 
evolutions. The plots in Fig. [TJ are therefore qualitatively 
valid for any R+aR~ n model (provided F > 0). It is also 
interesting to observe that the EOS is strongly varying 
with time just near z = 0: this should serve as a reminder 
that a simple parametrization of wde may often fail to 
describe interesting cosmologies. 

We can show now that an effective EOS w e s = 1/3 dur- 
ing the 0MDE is cosmologically unacceptable. In prin- 
ciple this should be shown case by case by a complete 
likelihood analysis of CMB and LSS data (see [16J for 
such an analysis for various coupled models), but this 
program is hardly feasible if we want to make general 
statements on f(R) theories. Instead we take a simpler 
but general approach. We calculate the angular size of 
the sound horizon 

9s -l^^w'k wr (14) 

where c 2 (z) = l/[3(l + 3ph/4p 7 )] is the adiabatic baryon- 
photon sound speed. According to the WMAP3y results 
[15j |. the currently measured value assuming a constant 
w is 9 S = 0.5946 ± 0.0021 deg. As radiation follows the 
same conservation law, the thermal history is the same 
as in usual cosmology so that Zdec is unchanged. It is 
easy to show that the integrand dz/H(z) is conformally 
invariant. For the model fi"2"Tl we integrate numerically 
the equations of motion in EF (including radiation) by 
changing initial conditions for x\, x%, x% via a trial and er- 
ror procedure until we obtain a present universe with a JF 
matter and radiation densities as observed in the WMAP 
data (we used 11™ = 0.3 and h = 0.7). Once we have 
the full background solution we evaluate 9 S with Zdcc ob- 
tained by solving the relation Zdcce K ^~^ ^^ = Zdcc- We 
have two competing effects: first, since the (j)MDE is more 
decelerated than in the usual case r(zdcc) will be system- 
atically smaller; second, the physical sound horizon dis- 
tance at decoupling is smaller than in usual cosmology 
partly because in our models a oc t 1 / 2 also between t eq 
and idee with z cq > Zdcc and mostly because H (z) before 
decoupling is much higher than in the standard case. As- 
suming H(z) = H (l + zf instead of (1 + zf' 2 , H(z dec ) 
becomes 30 times larger than in standard models. We 
find that the second effect is by far the dominating one, 
and as a consequence 9 S turns out to be an order of mag- 
nitude smaller than the observed value. In practice, we 
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Figure 2: The sound horizon angular distance 9 S as a function 
of the coupling /3 for n = 1 (thick line) and n = —2, 3, 10 
(dotted lines). The disk marks the value for /3 = 1/2. The 
grayed region shows the WMAP3y constraint at 4o\ 

find that 9 S can be approximated to a few percent using 
a standard cosmological model with an uncoupled dark 
energy component and a matter component with effec- 
tive equation of state w = 1/3. The typical values we 
find are near 9 S ~ 0.03 deg, i.e. more than ten times 
smaller than in a standard model. The periodic spacing 
A£ between the acoustic peaks in the CMB will be larger 
too by nearly the same factor. 

In Fig. [2] we plot the value of 9 S function of the 
coupling constant /3 (see eqs. [MS]) and for several n's. 

Clearly there is no way that small changes in fl^ or h 
or w c ff can cure this problem. Discrepancies are found 
as well for the age of the universe which turns out to be 
near 10-11 Gyr. As to be expected, we also find that 
the perturbations depart significantly from the standard 
case (as in pure exponential coupled models [IB]): for 
the matter density contrast on sub-horizon scales we find 
S ~ a 2 during 0MDE instead of the standard linear law. 

It is possible to generalize our results in several ways. 
First, one can show by direct substitution that the stan- 
dard matter era t 2 ^ 3 is a solution of (f2]) only for pure 
power laws (plus possibly a cosmological constant) R~ n 
with n = —1, — (7± v / 73)/12. In the last two cases, how- 
ever, the "matter" era occurs for Q, m = and is therefore 
unacceptable. The first case is clearly the pure Einstein 
case: this shows that a standard sequence of (exact) t 2 / 3 
expansion followed by acceleration can occur only for 



ACDM. In contrast, the c/>MDE generically exists as a 
saddle point. For R~ n with -1 < n < -3/4 the 0MDE 
is instead a stable point and the models are ruled out 
anyway. Still, this alone does not guarantee that the 
0MDE is always reached, regardless of the initial con- 
ditions, and a case-by-case numerical analysis is neces- 
sary. We explored extensively models like f(R) — exp(R) 
and log(i?) and always found </)MDE before accelera- 
tion. We also carried out a preliminary analysis of La- 
grangians for the models like / = R — c(i? /jy i?'" y ) _n , 
/ = R-c{R aMS R a ^ s y n and / = i?-c(i? GB )~" (where 
Rqb is a Gauss-Bonnet term) and did not find any ac- 
ceptable cosmological evolution. 

For the f(R) = R— fJ-f/R + H2R 2 models, two mech- 
anisms that could satisfy the local gravity constraints 
were suggested. One [6| achieves a short interaction 
range (or a large field mass) by adjusting ^1,^2 so that 
d 2 V(4>)/d(j) 2 vanishes today, when R = v3/if. Before this 
the R 2 term dominates and therefore we are back in one 
of our cases and the <^MDE takes place. Moreover, we 
find that the local minimum in the EF potential for such 
a class of Lagrangians does not lead to a late-time (effec- 
tive) cosmological constant. Another possibility to build 
a large mass is to take a very small H2 H , but in this case 
it is the 1/R term that dominates the cosmological evo- 
lution from the end of radiation, and again we are back 
in one of our cases. So even models that are designed to 
pass local gravity experiments fail our cosmological test. 
In summary, the main feature of our analysis is the mod- 
ification of the standard matter dominated epoch for the 
f(R) dark energy models investigated here. Hence these 
models are ruled out as viable cosmologies even if they 
are arranged to pass the Supernovae test and the local 
gravity constraints. We conjecture that our results apply 
to a much larger class of f{R) models; the precise con- 
ditions that determine the cosmological behavior will be 
published in future works. 
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